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Abstract Lame equation arises from deriving Laplace equation in ellipsoidal coor- 
dinates; in other words, it's called ellipsoidal harmonic equation. Lame functions are 
applicable to diverse areas such as boundary value problems in ellipsoidal geometry, 
chaotic Hamiltonian systems, the theory of Bose-Einstein condensates, etc. 

In this paper I will apply three term recurrence formulafll] to the power series 
expansion in closed forms of Lame function in the algebraic form(infinite series and 
polynomial) and its integral forms including all higher terms of A„'s. I will show 
how to transform power series expansion of Lame function to an integral formalism 
mathematically for cases of infinite series and polynomial. One interesting observa- 
tion resulting from the calculations is the fact that a 2^1 function recurs in each of 
sub-integral forms: the first sub-integral form contains zero term of A' n s, the second 
one contains one term of A„'s, the third one contains two terms of A,,'s, etc. Section 
6 contains additional examples of application in Lame function. 

This paper is 6th out of 10 in series "Special functions and three term recurrence 
formula (3TRF)". See section 7 for all the papers in the series. Previous paper in series 
deals with the power series expansion of Mathieu function and its integral formalism 
lfT4l . The next paper in the series describes the power series and integral forms of 
• . Lame equation in the Weierstrass's form and its asymptotic behaviors lfTBI . 
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1 Introduction 

The sphere is a geometrical perfect shape, the set of points which are all equidistant 
from its center (a fixed point) in three-dimensional space. In contrast, an ellipsoid 
is a imperfect one, a surface whose plane sections are all ellipses or circles; the set 
of points are not same distance from the center of the ellipsoid any more. As we all 
recognize, the nature is nonlinear and imperfect geometrically. For the purpose of 
simplification, we usually linearize those system in order to take a step to the future 
with a good numerical approximation. Actually, many geometrical spherical objects 
(earth, sun, black hole, etc) are not perfectly sphere in nature. The shape of those 
objects are closely better interpreted by an ellipsoid because of their rotations by 
themselves. For example, the ellipsoidal harmonics are represented in calculations 
of gravitational potential[9|. However spherical harmonic is preferred over the more 
mathematically complex ellipsoid harmonics (the coefficients in a power series ex- 
pansions of Lame equation have a relation between three different coefficients). 

In 1837, Gabriel Lame introduced second ordinary differential equation which 
has four regular singular points in the method of separation of variables applied to the 
Laplace equation in elliptic coordinates[2|. Various authors has called this equation 
as 'Lame equation' or 'ellipsoidal harmonic equation' ]5). 

Previously, there was no analytic solution in closed forms of Lame function j6][7] 
O . Using Frobenius method to obtain an analytic solution (represented either in the 
algebraic form or in the Weierstrass's form), the solution automatically comes out 
3 term recurrence relation(0|8]. In contrast, most of well-known special functions 
consist of two term recursion relation (Hypergeometric, Bessel, Legendre, Kummer 
functions, etc). 

In this paper I'll construct the power series expansion of Lame function in closed 
forms analytically and its integral forms with three-term recurrence formula [ 1 ). Lame 
equation is a second-order linear ordinary differential equation of the algebraic form[2| 



dx 2 2\x — a x — b x — cjdx 4(x — a)(x — b)(x — c) 

Lame's equation has four regular singular points: a, b, c and °°. Assume that its solu- 
tion is 




(1) 



y(z) = £ c„(x-a) 



n+X 



(2) 



n=0 



Plug © into ©■ 



c„+i =Ai c„+B n c„_i 



■n > 1 



(3) 



where, 



\(a(a + l)a-q) - {2a-b-c)(n + X) 2 
{a-b)(a-c){n + l+X){n + \+X) 
[a -(I -2(n+X))][a-2(n- 1+A)] 
2 2 (a -b)(a -c){n + l+X){n + \+l) 



(4b) 



(4a) 



c\ = A c 

We have two indicial roots which are Ai = and %i = \ 



(4c) 
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2 Power series 

2.1 Polynomial in which makes B n term terminated 

In Ref.[ 1 1, the general expression of power series of y(x) for polynomial of x which 
makes B„ term terminated is 

y(x) = £ c n x" +x = £ y m (x) = t c ? xM =yo(*) +yi(x) +y 2 (x) +y 3 (x) + ■■■ 

n=0 m=0 1=0 

f Po /ia-l \ 

=co\ i n^t+i U 2 

[ i =0 V'1=0 / 

Po ( io-1 ft (h-\ \ 1 

'0=0 t il=0 '2='0 \'3='0 / J 

°° r Po ( io-1 AT-1 / ft \ 

+e i e s A 2, n fi 2 !1+ in e a 2, 2 ,+a- n ^,+^+1) 

N = 2 l '0=0 I '1=° *=1 \'2t='2(<:-l) 'M+l=i2(*-l) / 

*. e ( 2 ff fi2i W+1+( ^ 1) > )}W jv+ 4 (5) 

'2JV-'2(«-l) \'2;v+i-' 2 («-i) / J J J 

For a polynomial, we need a condition, which is JTJ: 

B 2ft+(i+i)=0 where z',j3; = 0,1, 2,--- (6) 

In this paper Pochhammer symbol (x)„ is used to represent the rising factorial: (x) n — 
r(jc) • ^ n arJOVe ' A' i s an eigenvalue that makes B n term terminated at certain value 
of n. (|6]l makes each y, (x) where i = 0, 1,2, — as the polynomial in (fSJ. 



2.1.1 The case of a — 2 (2a,- + i + A) where i, a,- = 0, 1,2, • • • 

In d4a1i-(l4cb replace a by 2 (2a,- + / + A). In (O replace index /3, by a,. Take the new 
(I4al>-(l4cli. © and put them in (0 with replacing variable x by z. After the replacement 
process, the general expression of power series of y(z) for polynomial in which B„ 



4 



Yoon Seok Choun 



term is terminated is 



y(z) = L >^( z ) 

m=0 



and 



g (-ofr), (ao + ± + A) io ^ n 

i =0 (1 + 2")<o(4 + 2 )'0 

<^2 „/■„ i x 



(2a-fe-c)(/ + ^) 2 - a (a + ^)(ao + i + f) + ^ (-^(ap + I + A),- 

(*0 + 2 + §)(*> + ? + §) (l + T)<o(! + |).b 

(-«i),- 1 («i + | + A),-,(| + 4) i0 (| + ^) 

(-ai)ib («i + ! +*);o (1 + |k(f + 




' S (2 a -fc- c )(/ + |) 2 - a (a + |)(ao + | + |) + ^(-ao), (a + i + A)„ 
; =o ('0 + 3 + f X'o + l + f) (l + |)i (! + |)'b 

«* (2 fl -/7- c )(4+| + |) 2 - fl (a, + f + + f + i + |) + |r 



rii l 

(-a,), t (a, + fc + 1 + A)fr(l + 1 + (| + | + IV, ) 



2 



f (-«„),„ (a n + n + \ + A),„(l + § + §),„_, (| + § + §),„_, „ it ,, , 
.4_ i (- an ) (n _ 1 (a n +«+i+A) iB _ 1 (l + l + |) in (| + f + |) in f 



where 



Z = x — a 



T _ (fl-A)(a-c) ( 8 ) 
^ ~~ (a-6)(fl-c) 



a = 2(2(Xi + i + A) as j, a,- = 0, 1 , 2, ■ ■ ■ 
a, < a ; - only if i < j where i, j = 0, 1,2, 



(9) 



2.7.2 77ze case of (X — —2(2(Xi + i + A) — 1 where i, a, = 0, 1,2, • • • 

In d4ab-(l4cb replace a by — 2(2a,- + / + A) — 1. In (O replace index j3; by CU,. Take the 
new d4ab-d4ct. © and put them in (O with replacing variable x by z. Its solution is 

exactly equivalent to (Q. Put co= 1 as A=0 and cq = ( — 7 ] as A = A in 

\{a — b)(a — c) J A 

(IT). Then, we obtain two independent solutions of Lame equation. The solution is the 
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following ways. 



( \ (Ct \ 1 /a 1 \ 



— (x — a) 

' M = — 

/ o.i (a-b)(a-c) 



(a — b)(a — c) 



= g (-«o), (ao + l) iO T7 , n | g / (2 a -fc-c)/g- a ao(ao + I) + ^(-a ), (q + i), 



'0=0 




(4)10(^)10 i'o=0 

(-ai)( (ai + |)io(|)i 1 (f)i 1 

g (2a-fc-c)i§-aab(ab + ^) + ^ (-«o) io («o + ^), 

i =0 




(ztX'oWi'o 



('o + 5)('o + z) 



(4)('o(l)(0 



L 

l k='k-l 



(2a -b- c)(i k + |) 2 - a(a k 



§)(<%■ 



F 



(if 




, (-^k(^+fc+i), t (i + iV 1 (| + §k- 

' {-OkX-t (ak+k+Di^ (1 + 1),-,(| + |) ( - t 
' 4, (-«„)<„_, (a„ + » + \) in _ x (\ + f),„(| + §),„ 



77'" s/i" 



(10) 




1 /a \ 
2(2 +1+J ) 



— (x — a) 



y=0,l,2,-' (a-b)(a-c) 



f (-ao), («o + |), ^, n 
i'o=0 (4)% (l)'o 

(2 a -/7-c)(/o + |) 2 - a (ao + |)(ao + |) + ^ (-ao), (ao + |) io 

(*o + !)(*> + 2) (Ik(l)ib 
(-«i)i,(«i + |),-i(|)ib(|)ib ,-, 

S (2 a -fc-c)(/o + |) 2 - a (ao + i)(ao + i) + ^ (-ao) io (ao + |) io 

,o=o " ('0 + §)(*> + £) ~~ (!) io (l) i0 

a* 
I 



a * (2a-b-c){i k + k 1 + \) 2 -a{a k + k 1 + \){a k + k 1 + l 1 ) + f A 



x 



^=1 \ i k =i k -i {i k + j + l){i k + j + j) 

(-Ok)i k (ak + k + l)i k $ + §)/,_, (1 + I 
(-«*);*_, {a k +k+l) ik _, (| + §) it (l + \\ j 

f (-«»),„ (| + §),„_, (1 + §),„_, n ) 1 

\„=t_ 1 (-«nk-,(«n+« + lk-t(! + !k(l + fk J J 
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( [Tol l is the first kind of independent solution of Lame function for the polynomial as 
a — 2(2(Xj + j) or — 2(2ot/ + j) - 1 where j, 05/ = 0, 1,2, • • • . And (fTTT i is the second 
kind of independent solution of Lame function for the polynomial as a — 2(2a ; - + 
./') + 1 or -2(2a/ + ./' + 1 ) where 7, a,- = 0, 1 ,2, • • • . 



2.2 Infinite series 

In Ref.[ 1 1, the general expression of power series of y(x) for infinite series is 



y( x ) = E =yo{x)+y\{x)+y 2 {x)+y i {x) + --- 

m=0 

[i =0\<l=0 / '0=0 I il=0 i 2 =i \'3=i0 /J 

oc r oo r i -i w-i / oo i 2k -i \ 

+ E ^ E ^ A 2i Il^M + l J! E A 2, 2 ,+A: B 2i 2t+1 +(A:+1) 

^=2 L '0=0 I, l'l=0 k=l \'2/t='2(/i-l) '2t+i= 1 2(*-l) / 

x e f n ^ +1+ (, +1) )}l 2 - + ^} d2) 

l 2N=>2(N-\) \'2A'+1='2(/V-1) / J J J 

In ( fT2l replace a variable x by z and substitute (f4ab-d4cb into d!2l . Then, the general 
expression of power series of y(x) for infinite series is 



y{z) = E 3 '" 1 ^) =3 ; o(z)+yi(z)+3'2(z)+>'3(z) + --- 

m=0 

= coz a f £ (zi±iMi±i±ik Tj! -o + £ f (2fl-*-c)('o+_) 2 -£«(« + i)+ft 

[ I0 =0 (* + _)/o(!+ _)*o ( =0 1 0'o + 3 + §)(*0+4' + §) 

x + 2")'o(? + 4 + ?)'0 ^ J (~7 jl 2 jj 2")'' + 4 + 2")'! (2 + 2")'o(4 + 2")'0 „t| 

(1 + 2")i'o(4 + 2")('o i'i=!'o [ ( — 7 + 2 + 2") ! 'o(j + 4 + l) '0(2 + 2)11 (? + 27<1 

(2 a -fc- c )(/ + |) 2 -^a(q + l) + ^(-f + |), (f +i+|) i0 
('0 + I + §)(<<) + 5 + (1 + 1)4,(1 + 1)/ 

(2a - - c)(i k + \ + |) 2 - fra(a + 1) + £ 




^ & + § + _ + _)&+_ + ?+_) 



1 4 2! 2 ' 2 v 4 2! 2 ' 4 2! 2 /'t v 1 2! 2 ' 2 ^2 ' 4 ^ 2 ^-1 I 
f__.4-_4.Ay f24.l4.i4.Av fi 1 i 1 lv f- + - + -V I 

v 4 ' 2 ' 2 /'t-l W^2'4'2 V 1 ' 2 ' 2 > l k\2 ' 4 T 2 > l k ) 

{—« 4. _ J- Ay f« J- _ 4. I J. Ay fl J. _ 4. Ay f_4__4-Ay 

E : 4T 2 T 2 ;i„l 4 ^2 r 4 r 2/'^ 1T 2^ 2/'h-iV2 ' 4 2! 2/'«-l „ 
f_S4.54.Av fa4._4.I4.Av fl 4 « 4 AV f" 4 2 4 A V 



j~" (13) 
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( -1 \ k 

Put cq= 1 as A=0 and cq = ( -7 ) as A = A in fl3[ . Then, we obtain two 

\(a — b)[a — c) J 

independent solutions of Lame equation. The solution is the following ways. 

— (x — a) — (x — a) 2 



y(z) =LF[pi 



(a — b){a — c) ' (a — b)(a — < 
(-f)fc(f + ?)*,„«, , f f(2<i-&-c)ig-£a(a + l) + £ 





i =0 (4)i'o( 1 )'0 i =0 [ ('0+2)('0 + 4) 

x (~ "4 )'o( 7 + 4)10 ^ 1 (~7 j[ ijj] ^1 — llli CgjjoC?)^ 
(jJi'ovljiO i'i=!'o i ( — 7 "I" 2)10(7 4)k>\l)h 

(2a - fc - c)f§ - ^a(« + 1) + ft (-f ), (f + i), 
(io + + (!)*,(!)*, 
^ (2a-fr-c)(i t +§) 2 -fta(a + l) + ft 

4=4-1 ( I *+2 + 2)( l *+2 + ?) 

( _ 4 + 2)4-1 (? + 4 + 2)4-1 (1 + 2)4(4 + 2)4 j 

x ^ (~7 + 2)'" (j + 4 + fkiCj + 2)'n-l (4 + 2)'«-l j 
4=4-1 (~ 7 + f )'n-l ( 7 + 4 + i)4-l (1 + 3)^(4 + I)'" J 



(14) 



y(z) = LS hi = 



— (jc — a) — (x — a) 2 

— ; f] = — — 

(a — b)(a — c) (a — b)(a 





b-c){i + \f~%a{a + \)- 



{k + \){k 



( 4 + 4)4 (4 + 1)4 (3)40(2)* ^4 



+ 4)'o( 4 "t" 1)10(4)4 (2)4 
\2 _ ^n/^ _i_ 2) J- -£. 




)«o( 4 + 2)'0 



' + 2^4)4(4 + 2^2)4(4^ 2)4-1 (1 + 2)4-1 1 
( — 7 + 2 + l)4-l (f + 2 + 3)4-1 (4 + 2)4(1 + 5)4 J 

L (~f + § + |)4(7 jljjt ^3 it j)'"-' (j ^ j)'"-l n in 

(_ a,n ,ly fa.nilV (5 , »Wi i"V ' 



(15) 
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(TBI is the first kind of independent solution of Lame function for the infinite series. 
And ( TTSb is the second kind of independent solution of Lame function for the infinite 
series. 



3 Asymptotic behavior of the function y(z = x — a) and the boundary condition 
for x 

3.1 Infinite series 

Now let's test for convergence of infinite series of the analytic function y(z). As n 
goes to infinity, d4ab and d4bb are 

-(2a-b-c) -1 
limA n =A = v ( limB n =B = - — (16) 

ns>i [a — b)(a — c) «>i {a — b){a — c) 

Substitute (fT6b into (01. For n = 0, 1,2, ■ ■ ■ , it give 



CO 








Cl 


= AC 






C2 


= (A 2 4 


-2?)Cb 




C3 


-(A 3 4 


- 2AB)C Q 


c 4 


= (A 4 H 


- 3A 2 B + 


B 2 )C 


C5 


= (A 5 H 


-4A 3 B + 


3AB 2 )C 


C6 


= (A 6 H 


-5A 4 B + 


6A 2 B 2 +B 3 )C Q 


C'7 


= (A 7 H 


-6A 5 B + 


10A 3 B 2 +4AB 3 )C 


C8 


= (A 8 H 


-1A 6 B + 


15A 4 B 2 + 10A 2 B 3 +B 4 )C 



(17) 



The sequences c„ consists of combinations A and B in ST% . First observe the term 
inside parentheses of sequence c„ which does not include any A„'s: c„ with even 
index (cq,C2, 04,- ■■). 
(1) Zero term of A's 

co 

c 2 = Bc 
C4 = B 2 cq 

c 6 =B 3 cq (18) 

eg = B 4 c 

c lQ =B 5 c 
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When a function y(z), analytic at z=0, is expanded in a power series z=0, we write 

y(z) = £ y m (z) (19) 

m=0 

where 

y»(z) = (20) 

n=0 

Put(fT8l> in ([20]) putting m = 0. 

yo(z)=c £(Bz 2 )" (21) 

Observe the terms inside parentheses of sequence c„ which include one term of A n 's 
in dTTb : c„ with odd index (ci, C3, cs,- • • ). 
(2) One term of A's 



C\ 


= Ac 


C3 


= 2ABc 


C5 


= 3AB 2 c 


Cl 


= 4AB 3 c 


C9 


= 5AB 4 c 



(22) 



Put (f22| in (|2(2l putting m= 1. 

yi(z)=coAz£^±H(Bz 2 ) n (23) 

Observe the terms inside parentheses of sequence c„ which include two terms of A„'s 
in dTTb : c„ with even index (C2, C4, eg, - • • ). 
(3) Two terms of A's 

c 2 = A 2 c Q 
c 4 = 3A 2 Bc 
c 6 = 6A 2 B 2 c 

c & = 10A 2 B 3 co W 
cio = l5A 2 B 4 c 



Put (|23]) in <[20j» putting m = 2. 



y 2 (z)= C o(Az) 2 f: (w + 1 )i W + 2) ( fe 2 )" (25) 

n=0 L - 
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By repeating this process for all higher terms of A's, I obtain every y m (x) terms where 
m > 3. Substitute (ETT i. (l23l . (f25l and including all y„,(x) terms where ra > 3 into (1191 1. 

y(z) = £c nZ " +A =y (z)+3'iW+3'2W+3'3(z) + --- (26) 

n=0 

= E E (w |^ )! fy" = — -i where c = l,x = Bz 2 and y = Az 

Substitute (O in ( 1261 

Hri y(z) = : 1 7 a-.? ] g^E^ T (27) 

1 ~r \(a-b){a-c) T (a-b)(a-c) J 

d27T i is geometric series. The condition of convergence of it is 

(x — a) 2 (2a — b — c)(x — a) 



(a — b)(a — c) (a — b)(a — c) 



< 1 (28) 



According to (1281 1. we have a boundary condition of x for the infinite series of Lame 
function in which is following the way. 

(A) As a=b or a=c 

no solution (29a) 

(B) As < ( fl -^y^) 2 < (a-b)(a-c) 



(b + c) f/ b + c\2 " ~ " (b + c) ft b + c\2 " ~~ " 

21 V V ~^r) +i > a - b ^ a -^ <x< — 1 -^ + \j\ a -—) +(a-b)(a-c) 

(C) As (a-b)(a-c) = ( fl -^y^) 2 >° 



(29b) 



{b + c) 



V2( fl -fc)(fl-c)<i<^ and ^-^<x<^-^ + ^2(a-b)(a-c) 



2 v v 2 2 2 

(29c) 

b + c\ 2 



(D) As < (a - b) (a - c) < la — 



(fl-^) 2 + {a-b)(a-c)<x< ( ^-J[a-^ 2 -(a-b)(a-c) 



i± T 1 + ]]( a - t T) 2 ' -{a-b){a-c)<x<^ + s J(a- b -^ 2 ' + (a-b)(a~c) 

(29d) 

(E) As < (a— <-{a-b)(a-c) 

{b + c) \~( b + c\2 " ~~ " (b + c) [} b + c\2 " ~~ " 
— 2 VV 2/ -\ a ~ b K a - c ) <x< — 2 h Vv 2/ -(a-oj(fl-c) 



(29e) 
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(F) As -{a-b){a-c)= (a-^-^ >0 

(b + c) r— —7 r (b + c) , (b + c) (b + c) r^rr- 

y - - y/-2(a-b)(a-c)<x< - and v - ' <x< y - - + ^-2(a- 

(29f) 

(G) AsO < -(a-b){a-c) < (a-^^j 

it^-^(a-^y-(a-b)(a-c)<x<^-^(a-^ 2 + (a-b)(a-c) 

[± l £l + \J( a -' 2 ¥) 2 ' + {a-b){a-c)<x<^ + ^(a- b -^f -(a-b)(a-c) 

(29g) 



3.2 The case of 2a — b — c^$> 1 or 2a — Z? — c <§; — 1 

Let assume that 2a — ^ — c is much greater than 1 or is much less than — 1 . Then B„ 
terms are negligible. (|3) is approximately 

c„ + i«A„c„ (30a) 

And, 

— (2a — b — c) 

hmA„ rs A = £ ci = Aqc q « Ac (30b) 

n»l (a — b)[a — c) 

Substitute ( I30bb into ( I30al >. For n = 0, 1 , 2, • • • , it give 

co 

c\ =AC 
c 2 =A 2 C 

C3 =A 3 C (3D 
c 4 =A 4 C 



When a function y(z), analytic at z=0, is expanded in a power series z=0 by using 
(l3Tl l. we write 

limy(z) =cqY (Az)" = — !— where |A Z | < 1 and c Q = 1 (32) 
n>l .r^, 1 - Az 



Substitute OObl into d32l) 



limy(z) = 7- -T- where a^b and a ^ c (33) 

«>1 1 1 ( (2a-b-c)(x-a) \ 

L + { (a-b)(a-c) ) 
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d3~3l is binomial series. The condition of convergence of x is 

(2a — b — c)(x — a) 



(a — b)(a — c) 



< 1 



(34) 



3.3 The case of 2a — b — c rs 

Let assume that 2a —b — c is approximately close to 0. But a ^ in (|5j- Then A„ 
terms are negligible. (|3) is approximately 



And, 



c n+ \ « B„c„_! where « > 1 



limB,, k,B = — 

n>l (a — b)(a — c) 



We can classify c„ as to even and odd terms from plugging (I35bt into 



(35a) 



(35b) 



CO 




C] 




C2 = 


Bc 


C3 = 


Bc\ 


C 4 = 


B 2 c Q 


C5 = 


B 2 Cl 


C6 = 


B 3 c Q 


Cl = 


B 3 Cl 


C8 = 


B 4 c 


C 9 = 


B A Cy 


ClO = 


--B 5 c 


cn = 


= B 5 Cl 



(36) 



When a function y(z), analytic at z=0, is expanded in a power series z=0 by using 
we write 



linWz) =c t (Bz 2 )" + C iz E ( Bz 



2\" . 



1 



l-Bz 2 ' Cl 'l-Bz 2 



n=0 



where|Bz 2 | < 1 



Substitute d35bt into ( l37b . And for simplicity, let say eg = c\ = 1 into it. 

1 + (x — a) 



limy(z) 



1 i ( (x-a) 2 
1 T \{a-b)(a-c) 



where a ^ b and a^c 



(37) 



(38) 



The condition of convergence of x is 



(x — a) 2 
(a — b)(a — c) 



< 1 



(39) 
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4 Integral Formalism 

4. 1 Polynomial in which makes B„ term terminated 

Now let's investigate the integral formalism for the polynomial case of B n term ter- 
minated at certain eigenvalue. There is a generalized hypergeometric function which 
is: In this paper Pochhammer symbol (x) n is used to represent the rising factorial: 

r(.v+n) 



h= f (-«/),-,(»/ + 1 + 1 + 1 +h)u-i (|+|+ jVi ,-, (40) 

- B(/,_ 1 + ^-I + |,./ + l)B(/ / _ 1 + i + |,./ + l)(/ / _ 1 -a / ) / (a / + /,_ 1 + Z + ^ + A) / . 

= TJ ) - - : : \ ^ : : - 1 

By using integral form of beta function, 

BL x + L-\ + \,i+l)= fd tl ^-^{l- tl y (41a) 



+ ~ + + = frfM/M?- 1 ^ 1+ * (!-»,)■' (41b) 



Substitute (14 lab and (!41bb into (flOt , And divide (z'/_i + 1 - | + §)(z'/-i + | + §) into 



' A-i (-«/)/z-i(«; + ^ + | + ^)iw(l + 3 + !)<;(! + 2 + I)''/ 



x £ ft-i -o>)j(ii + + _ _ j 



The integral form of hypergeometric function is 



where Re(y- a) > 
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replaced a, j3, y and z by z'/_i — 0C/, 2/_i + Z + 0C/ + | + A, 1 and 77 (1 — f/)(l — «/) in 



2niJ~" l Vi y ~ " ,/v " ^l-T]v/(l-f/)(l-M/)J 

(44) 



y'=o 

^v / -(l-J7v/(l-f/)(l-«,))- ( ' + ' +X) ' 



1 



\iX-±.){\- m {\- tl ){\-ui)) 

Substitute (0U) into (l42l >. 

1 

('/-l + 5 ~ 4" + |)0'z-l + 5 + 7) 



x 3 (-«/),-, (g + / + 1 + A)„(l + 1 + |v, (| + | + IV, 

(-«/)«'« (<* + ! + k + A )'/-i (! + 1 + + 2 + §)>i 

/ 3+7 / 1+7 — i dv,-(l-i7v,(l -„,))-('+?+*) 

Jo Jo 27n J v ; 



1 1 \ / \ ii 1 

1 - - \ / f /M ,v, 77 



(45) 



^l-T7v/(l-*z)(l-wz) y V^- 1 ) l-i?v/(l-^)(l-M Z ), 

Substitute d45l > into (0 where / = 1,2,3, and the integral formalism of Lame 
function for polynomial in which B„ term terminated is 

y(z) = E - Vm ( z ) 

m=0 



L '0=0 K L + 2M4 + 2>H> 
n-1 / r \ 



J0 



-A 



— — <£ dv n - k — (1 - V„_jt+i „V„-jfc(l -*„-*) (1 -Un-k)) { " k+1+X) 

2m J v„- k 

V n -k 1 — ^n-t+l, n V„-t(l — f„_fc)(l — U n -k) 



n-ivi / n-k.n 

x £ (-^ + l+A),o V ,l „| (46) 

i =0 (1 + 2)'o(4+ 2)'0 J J 



Lame equation in the algebraic form 



15 



where 



1 '»•'; : /'•»//»» 

Vy - <j (V, - 1) l-WljV ; (l-f/)(l-««) 

77 only if z > /' 



(47) 



Put cq= 1 as A=0 and cq = 



1 



, as A = i in d46l i. Also apply d43b into 

{a — b)(a — c)J 



it. Then I obtain two independent solutions of Lame equation. The solution is the 
following ways. 



l ra, 



y(z) = LF aj \aj = 

— (x— a) 2 



1 /a 1 



/ lor 1 1 h ; 

2V2 / 2V2 2 



;m = 



-(x — a) 



y=0,l,2,- (a-b)(a — c) 



(a — b)(a — c) 



1 3 



n— 1 / /-l 



2 Fi ( -oq, ao + -; -; J] ) + £ <j ]~[ ( £ df„_* tf_ k J> £ m^ a . 



= i /,-=o 



(fdv n -k-^— (l - Wn-ft+i „V n _fc(l -*„-jfc) (1 -M»-Jfc)) ( " /:+?) 
27H J V„_& 

(v„-jt-l) 1 



V«-A 1 - W„_jt + i „V„_i(l — f„_^)(l - H„_jt) 



c I (2a-* - c)Kit fc_1) (^-M^ n _ M ) 2 V^- 
-a f Ofe_jt-i + -(n— 1) J (a M + \(n-k-\y\ + 



x 2 ^i ( -ao,«o + i;^; Vi,„ ) ^ju" 



(48) 
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y(z)=LSa J [a j = -(----j)or--(- + l + j 

— (x — a) 2 \ 



— (x — a) 



;=0,1,2,- (a — b)(a — c) 



(a — b)(a — c) 



3 5 



= r?M 2*i ( -ao>ao+-;-;T? 



oc r «-i / ^ 

in 

«=1 I <r=0 



,*(»-*-2) 



du^h it 



n-k u n -k 



r("-*-4) 



— « <iv„_* (1 -V„-t + l,„V„-l(l -f„-it)(l -Un-k)) 

2m J v„_ k 



(n-k+i) 



{Vn-k~ 1) 



1 



1 - V„_ J t + i i „V„_ jt (l - f„_ t ) (1 — , 

2 - -i^B-t-i) 



(2a-Z>-c)Kir 4) 





1 




( <%i-*-l 






+ 2 





2 4 



(49) 



(l48l > is the integral formalism of the first kind of independent solution of Lame func- 
tion for the polynomial as a = 2(2(Xj+ j) or — 2(2ofy + j) — 1 where y, 0Cy = 0, 1 ,2, • • • . 
And, (l49l is the integral formalism of the second kind of independent solution of 
Lame function for the polynomial as a — 2(2aj + l or — 2(2a 7 - + j+l) where 
j,(Xj =0,1,2, .... 



4.2 Infinite series 



For infinite series, replace the finite summation with an interval [0, cto] by infinite 
summation with an interval [Q,°°] in (l46*l l. Also, replaced Oq, 0C„_a and tt n -k-\ by 
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^(!-4Kf-" + *- A ) -di(f-n+*+l-A)into®. 



y(z) = Ly«(z) 

n=0 

- ^7* J Y + |)'o(f + I + §)'0 fp 

1.10=0 l 1 + 2>'o\A + im 



n=l (. i=0 

X 



70 



— i dv n - k (1 - V n _ fe+ i^v„_jfc(l - ?„_*)( 1 - 

27H J v n _ k 

[Vn- k -l) 1 x|(f-« + ^) 



Vn-k 1 ~ ^n-i+ 1 ,n V n -k ( 1 ~ t n -k) ( 1 ~ , 

X ^ (2a - i - c) w \ w n - k . n d^. n W l_ k „ 



^ r <50) 



Put co= 1 as A=0 and co = (7 t-, r ] as A = A in eq(118). Also apply 

\{a-b)(a-c)J 

eq(99) into it. Then I obtain two independent solutions of Lame equation. The solu- 
tion is the following ways. 



-6c — a) —(x — a) 2 



(a — b){a — c) (a — b)(a — c) 

2fl hl'I + Z ; Z ; 1 +1 Ul /„ dt n-k Ck 2 / du «-k uf_ k 



4'4 44 / ^ iA \/n " _ * /n 

t t h- t / „ = i ^yt = o V 70 1/0 

;i ldv n ^k-^— (1 - V„_ fc+ i „v„_ fc (l -*„_*) (1 — ( " * +?) 

V n -k 1 — ^n-it+l,nV„-fc(l — f„_fc)(l - M„-jt), 



x«j (2a-fo-c)¥ 7 „4„ ^V n _ M dy n M J w „ 2 _ M - ^ «(« + 1) + ^4 
a a 13 
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/ — (x — a) —(x — a) 2 \ 

y(z) =US[n = \ > ;r,= - \ ' 

\ (a — b)[a — c) [a — b)(a — c) J 



a la 15 \ " [ / f 1 , A 



I? 3 <J 2 f l(-4+4'4+2 ; 4 ;r J 111 /„ dt n-kC-k L ""» • 



n=\ I *=0 



(n-k-2) f 1 5 (»-*-* 



JO 



X— ldv n - k -^— (l - V n _jt+l,nV n _ k (l -*„_*) (1 -«„_*)) ( " /:+t) 

(v„_,-l) 1 Nfc*-" ] 



Vn-i 1 — W'„_ J t + i,„V„_it(l — *„_*) (1 — Mn-fc), 

xi(2a~b-c)w n \ n ^ w n - M <V„-,„J w 2 n _ Kn - ^a(a + 1) + ^ > I 

**(-f4f4i^)H 

( l5TT l is the integral formalism of the first kind of independent solution of Lame func- 
tion for the infinite series. And, (l52l is the integral formalism of the second kind of 
independent solution of Lame function for the infinite series. 



5 Application to Laplace equation in ellipsoidal coordinates and nonlinear 
evolution equations 

A great many authors have worked on applications to boundary value problems for 
the Laplace equation in elliptic coordinates. We can apply the power series expan- 
sion in closed forms of Lame function and its integral forms into many mathemati- 
cal physics areas. For example, Lame equation can be employed to solve boundary- 
value problems for Laplace equation in elliptical cones. In "Occurrence of periodic 
Lame functions at bifurcations in chaotic Hamiltonian systems"|f3], the authors shows 
that Lame equation in the Weierstrass's form occurs at bifurcations in chaotic two- 
dimensional Hamiltonian systems with mixed phase space, (see (6), (7) in Ref.[3]); 
Lame equation in the algebraic form can be transforms to in the Weierstrass's form 
by changing an independent variable. Then we can describe Lame function more an- 
alytically in power series expansions in closed forms and its integral forms. In "Lame 
Function and Its Application to Some Nonlinear Evolution Equations"|5|, the Lame 
equation is applied to solve nonlinear (1 + l)-dimensional, (1 +2)- dimensional and 
coupled evolution equations, (see (1), (3), (5), (6) in Ref.0) 



6 Conclusion 

The power series expansion of Lame equation (ellipsoidal harmonics equation) has a 
relation between three different coefficients as we see ((3}- Because of its three term 
recurrence, there is a mathematical difficulty that Lame functions are described in the 
power series expansion in closed forms and its integral forms. 



Lame equation in the algebraic form 



19 



In this paper I derive the power series expansion in closed forms for the infinite 
series and the polynomial of Lame function analytically by applying three term re- 
currence formula[ l \. Even if coefficient a can have two different values in which are 
2 (2a, + i + A) or — 2 (2a,- + i + A) — 1 to make B„ term terminated at certain eigen- 
values, its solutions are completely equivalent to each other as we see Q. 

Also as we see representations in the form of integrals in Lame functions, a jF\ 
function recurs in each of sub-integral forms: the first sub-integral form contains zero 
term of A' n s, the second one contains one term of A„'s, the third one contains two 
terms of A„'s, etc. Then we can transform Lame functions to all other well-know 
special functions such as Hypergeometric, Bessel, Legendre, Kummer functions, etc. 

Since we get the integral forms of power series expansions in Lame function, 
we are able to obtain generating functions of it. The generating functions are really 
helpful in order to derive orthogonal relations, recursion relations and expectation 
values of physical quantities. 

Most of well-known papers with boundary value problems in ellipsoidal geom- 
etry have been published with using the Weierstrass's form of it. So we need know 
what the power series expansion in closed forms of Lame function in the Weier- 
strass's form and its integral forms. In lfl6l . I am deriving the power series expansion 
in closed forms of Lame function in the Weierstrass's form and its integral forms 
analytically. In ItTTl . I am investigating the generating functions of Lame equation in 
the Weierstrass's form. 



7 Series "Special functions and three term recurrence formula (3TRF)" 

This paper is 6th out of 10. 

1 . "Approximative solution of the spin free Hamiltonian involving only scalar po- 
tential for the q — q system" ifTUl - In order to solve the spin-free Hamiltonian with 
light quark masses we are led to develop a totally new kind of special function the- 
ory in mathematics that generalize all existing theories of confluent hypergeometric 
types. We call it the Grand Confluent Hypergeometric Function. Our new solution 
produces previously unknown extra hidden quantum numbers relevant for descrip- 
tion of supersymmetry and for generating new mass formulas. 

2. "Generalization of the three-term recurrence formula and its applications" ifTTI 

- Generalize three term recurrence formula in linear differential equation. Obtain the 
exact solution of the three term recurrence for polynomials and infinite series. 

3. "The analytic solution for the power series expansion of Heun function" ITLZl 

- Apply three term recurrence formula to the power series exp ansion in closed forms 
of Heun function (infinite series and polynomials) including all higher terms of A„s. 

4. "Asymptotic behavior of Heun function and its integral formalism", lfl3l - Ap- 
ply three term recurrence formula, derive the integral formalism, and analyze the 
asymptotic behavior of Heun function (including all higher terms of A„s). 
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5. "The power series expansion of Mathieu function and its integral formalism", 
|[T4l - Apply three term recurrence formula, analyze the power series expansion of 
Mathieu function and its integral forms. 

6. "Lame equation in the algebraic form" |[T5l - Applying three term recurrence 
formula, analyze the power series expansion of Lame function in the algebraic form 
and its integral forms. 

7. "Power series and integral forms of Lame equation in the Weierstrass's form 
and its asymptotic behaviors" [ 16 1 - Applying three term recurrence formula, derive 
the power series expansion of Lame function in the Weierstrass's form and its integral 
forms. 

8. "The generating functions of Lame equation in the Weierstrass's form" ifTTI - 
Derive the generating functions of Lame function in the Weierstrass's form (including 
all higher terms of A„'s). Apply integral forms of Lame functions in the Weierstrass's 
form. 

9. "Analytic solution for grand confluent hypergeometric function" Ifl8l - Apply 
three term recurrence formula, and formulate the exact analytic solution of grand 
confluent hypergeometric function (including all higher terms of A„'s). Replacing ji 
and E(0 by 1 and — q, transforms the grand confluent hypergeometric function into 
Biconfluent Heun function. 

10. "The integral formalism and the generating function of grand confluent hy- 
pergeometric function" lfl9l - Apply three term recurrence formula, and construct 
an integral formalism and a generating function of grand confluent hypergeometric 
function (including all higher terms of A„'s). 
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